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Abstract 

We report on the discovery of a transient X-ray pulsar, AX J1841. 0—0536, serendipitously found in 
the Scutum arm region with the ASCA in two separate observations. The X-ray flux is very faint at the 
beginning, but exhibits two flares in the second observation. The flare flux increases by a factor 10 within 
only '^1 hr. Coherent pulsations with a period of 4.7394±0.0008 s were detected in the brightest flare 
phase. The X-ray spectra in the quiescent and flare phases were fitted with an absorbed power-law model 
with a photon index ~ 1 plus a narrow Gaussian line at the center energy of 6.4 keV. The interstellar 
column density of ~ 3 x 10^^ cm~^ may indicate that AX J1841. 0—0536 is located at a tangential point 
of the Scutum arm at 10 kpc distance. The coherent pulsations, large flux variability and the spectral 
shape suggest that AX J1841. 0—0536 is a Be/X-ray binary pulsar. 
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1. Introduction 

The Scutum region near I 30° is in the tan- 
gential line of sight of a galactic arm (the Scutum 
arm), and hence shows intensity excesses in 
the infrared, radio and 7-ray regions (e.g., |^00 
[cite] cite. hayakawaHayakawa et al.(1977)) as well as 
diffuse X-rays (^0 [cite]cite.yamauchil993Yamauchi, 
Koyama(1993); ^0 [citeJcite.sugizakiSugizaki et 
al.(2001)). This region is also rich in X-ray binary 
sources; many transient X-ray sources have also been dis- 
covered (e.g., ^0 [cite] cite. koyamaKoyama et al.(1990); 
|00 [cite]cite.yamauchil995Yamauchi et al.(1995); gOO 
[cite] cite. teradaTerada et al.(1999)). From the spec- 
tral and flux variability, these X-ray transients are 
suspected to be either neutron-star binaries with high 
mass or low-mass stars, or black-hole binaries (|^00 
[cite] cite. tanakaTanaka, Shibazaki(1996)). Coherent 
X-ray pulsations have been found from some of them, 
suggesting Be/X-ray binaries. 

Among them, an unusual source is AX 1845.0—0433, 
which exhibited violent flares with a very rapid rise 
time of < 10^ s (^0 [cite]cite.yamauchil995Yamauchi et 
al.(1995)). Apart from the unusual flares, other charac- 
teristics are similar to those of high-mass X-ray binaries. 
In fact, ^Jcite.coeCoe et al.(1996) ([citeJcite.coeCoe et 
al.(1996)) reported an 09.5 star as the optical counterpart 
of AX 1845.0—0433. Nevertheless, no coherent pulsation 
has been found. 

In the ASCA survey and pointing observations on the 



Scutum arm region, we found another unusual source, 
AX J1841. 0—0536, which also exhibited violent and rapid 
flares (^00 [cite]cite.bambal999Bamba, Koyama (1999)). 
This paper reports on a detailed analysis of this new tran- 
sient X-ray source, particularly on the detection of coher- 
ent pulsations and phase resolved spectral analyses. 

2. Observations and Data Reduction 

The first ASCA observation of a region around {l,h) — 
(26.°9, -0.°1) was made on 1994 April 12, or MJD 
49454.677-49454.701, as a part of the Scutum arm survey 
project (hereinafter, obs. 1). A longer exposure follow-up 
observation was made on 1999 October 3-4, or on MJD 
51454.252-51455.463 (obs. 2). ASCA carried four XRTs 
(X-Ray Telescopes, ^0 [cite]cite.serlemitsosSerlemitsos et 
al.(1995)) with two GISs (Gas Imaging Spectrometers, |^00 
[citeJcite.ohashiOhashi et al.(1996)) and two SISs (Solid- 
state Imaging Spectrometers, |^00 [citeJcite.burkeBurke et 
al.(1994)) on the focal planes. Since our target was out of 
the SIS field in both observations, we do not refer to the 
SIS in this paper. The GISs were operated in the nominal 
PH mode with a time resolution of 62.5 ms (High bit-rate) 
or 0.5 s (Medium bit-rate). 

We rejected the CIS data obtained in the South Atlantic 
Anomaly, in low cut-off rigidity regions (< 6 GV), or when 
the target's elevation angle was low (< 5°). Particle events 
were removed by the rise-time discrimination method (|^00 
[citeJcite.ohashiOhashi et al.(1996)). After these screen- 
ings, the total available exposure times of obs. 1 and 2 
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were ~ 2 ks and ~ 36 ks, respectively. To increase the 
statistics, the data of the two detectors, GIS-2 and GIS-3, 
were combined in a following study. 

3. Results 

3.1. X-Ray Image 



source is found in both obs. 1 and obs. 
a point-like source image located 



A transient 
2. Figure shows 
near the galactic plane. The source position is deter- 
mined by the method of |c]citc.gotthclf2000Gotthclf et 
al.(2000) ([cite]cite.gotthelf2000Gotthelf et al.(2000)) and 
is (R.A., Dec.)j2ooo = (18'^41'"50?2, -05°36'00") with an 
uncertainty of ~ 1'. Therefore, we designate this source 
as AX J1841. 0—0536 (this source was at first named 
AX J1841. 0—0535.8 based on less-accurate positioning; 
|00 [cite]cite.bambal999Bamba, Koyama (1999)). No X- 
ray source has been found within the error region in the 
SIMBAD database.^ However, this source is m an error 
region (0.°4 x 4°) of a previous X-ray transient (the Ginga 
source No. 2 in 000 [cite] cite. koyamaKoyamaet al.(1990)). 




Fig. 1. X-ray image around AX J1849. 8-0530 in the 
2.0-7.0 keV band combined with both observations, with the 
galactic coordinates. G26.6— 0.1 is an unidentified hard dif- 
fuse source (A. Bamba et al. unpublished). The source and 
background regions for the spectral and timing analyses are 
also indicated by the solid and dotted circles, respectively. 



3.2. Light Curve and Timing Analyses 

The light curves of AX J1841. 0-0536 in both observa- 
tions are shown in figure ^ The X-ray flux is highly vari- 
able with two flares of the on-set epochs at MJD 51454.7 
and MJD 51455.3 with an interval of only 0.6 d. As 
can be seen in figure |^, the flares are found in the hard 
bands, while the soft band fluxes are essentially constant. 
A remarkable fact is that the on-set time of both the flares 
are extremely rapid; the fluxes increased by one order of 
magnitude within 1 hr. 
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Fig. 2. X-ray light curves ofAXJ1841. 0-0536 observed with 
GIS-2 and GIS-3 (bin time = 1600 s) in obs. 1 and obs. 2. 
The upper and lower panels are those in the 0.7—2.0 keV 
and 2.0-10.0 keV bands, respectively. For the spectral anal- 
yses (see text), obs. 2 was divided into three epochs (obs. 2a, 
obs. 2b, and obs. 2c), as indicated in the lower panel. 



To study the time evolution, we divided obs. 2 
into three epochs: obs. 2a (MJD 51454.252-51454.700), 
obs. 2b (MJD 51454.700-51455.300), and obs. 2c (MJD 
51455.300-51455.463), as shown in figure |. We then 
performed detailed timing analyses of the brightest flare 
episode (obs. 2c). We extracted ^ 1000 photons, including 
^ 90 background photons, in the 1.9-4.9 keV band from 
a 3' radius circle around the source (the solid circle in fig- 
ure 1), corrected the photon arrival times to those at the 
barycenter and applied an FFT (Fast Fourier Transform) 
algorithm. The resulting power-density spectrum, figure ^ 
(left), shows a peak at ~ 0.211 Hz. With 13471 frequen- 
cies examined, the chance probability of producing such a 
large peak is ^ 9.8 x 10~^, and hence the 0.211 Hz oscil- 
lation is highly significant. We next performed an epoch- 
folding search and determined a more accurate pulse pe- 
riod of P = 4.7394 ± 0.0008 s. Figure | (right) shows the 
folded pulse profile in the soft (1.9-4.9 keV) and hard 
(4.9-10.0 keV) bands. 

The pulse shapes are sinusoidal in both energy bands, 
but the background-subtracted pulse fraction is larger 
in the soft (28%) band than in the hard (8%) band. 
Although same analyses, FFT and epoch folding, were 
performed for the other data sets, no significant power 
was found with the 3a upper limit of the background- 
subtracted pulse fraction in the soft band of about 45% 
and 12% for obs. 2a and 2b, respectively. 

3.3. Spectral Analyses 

We made time-sliced X-ray spectra using photons from 
the same region of the timing analysis (in the 3' radius 
circle around the source; see figure 1). The contribu- 
tion of the galactic ridge emission is not negligible for 
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Fig. 3. Upper panel: Power-density spectrum in obs. 2c. 
A significant peak is detected around ~ 0.211 Hz. Lower 
panel; Pulse profiles in the 1.9-4.9 keV (upper panel) and 
the 4.9—10.0 keV (lower panel) bands in obs. 2c. The vertical 
axes indicate the normalized intensity per bin. The horizon- 
tal axis indicates the pulsation phase and the epoch of the 
beginning of the phase corresponds to MJD 51455.300. The 
background levels are indicated by the broken linos. 

AX J1841. 0—0536 located closely on the galactic plane 
at {l,b) = (26.°7, 0.°1). Contamination from an unidenti- 
fied hard diffuse source G26.6— 0.1 (A. Bamba et al. un- 
published) and vignetting may not be ignored. We thus 
made the background spectrum by adding data from two 
source-free regions, indicated by the dottcd-circlcs in fig- 
ure |l|. These two regions were selected to compensate 
for the latitude-dependent galactic ridge emission (|^00 
[cite] cite. kanedaKaneda et al.(1997)) and contamination 
of G26.6— 0.1. The vignetting is also the same as the 
source region. 

The background-subtracted spectrum in the brightest 
flare phase (the second flare, or obs. 2c) was at first fit- 
ted with a power-law function. The absorption column 
was estimated from the cross sections of Morrison and 
McCammon (1983) with solar abundances. Although this 
model is statistically acceptable, we see a systematic data 
excess at around 6-7 keV. Adding a narrow Gaussian line 
to the model reduced the from 51.2/52 to 43.2/50 
(x^/degree of freedom); thus, the presence of the fine is 
significant. The center energy of the line is 6.40±0.12 keV 
(hereinafter, the errors are of 90% confidence level, unless 
otherwise mentioned), implying that the line is due to 
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neutral or low-ionized irons. 

An absorbed thin-thermal plasma model fitting, on the 
other hand, requires unrcalistically high temperature of > 
80 keV, and is hence unlikely. 

We then fit the background-subtracted spectra of the 
other epochs with a power-law plus a Gaussian model, 
fixing the center energy of the line to be 6.40 keV (the 
best-fit value in obs. 2c). These fittings are acceptable 
for all epochs. The spectra and the best-fit models for 
the four separate epochs are shown in figure ^, while the 
best-fit parameters are listed in table 1. 

4. Discussion 

AX J1841. 0—0536 is in the error region 
(/ = 26.°6 ± 0.°2,6 = -0.°5 ± 2.°0) of a previous X- 
ray transient No. 2 in ^Jcite.koyamaKoyama et 
al.(1990) ([cite] cite. koyamaKoyama et al.(1990)). 
Although the position uncertainty is large, the X- 
ray features, power-law index of 1.5±0.3, flux of 
2 X 10"" erg cm'^ s'^ and log (Nh) of 22. 7± 0.3 are 
similar to those of AX J1841. 0—0536 in obs. 1 and 
obs. 2b. Therefore these two transients are likely to 
be the same source. From the flux history, we suspect 
that the quiescent level of AX J1841. 0—0536 is on the 
order of 2 x 10"^^ erg cm~^ s~^, and that flares may 
occur very frequently. More bright flares of a level of 
10"^'' erg cm^^ s~^ are, however not frequent. 

Although the limited statistics allows a constant-value 
hypothesis of iVn, we see a hint that Nn increase with 
increasing flux, which suggests that the circum-stellar ab- 
sorption becomes larger in the fiare, due to the larger gas 
accretion. Therefore, the interstellar absorption should be 
regarded as the minimum value of ~ 3.2t_2'^ x 10^^ cm s~^ 
(obs. 1 and obs. 2a). The large error of A^h does not allow 
us to reliably estimate the source distance. Nevertheless, 
we simply estimate the distance to AX J1841. 0—0536 to 
be 1-10 kpc, assuming the average density in the galactic 
plane to be 1 H cm^'^. This value is consistent with that 
the transient source is located in the tangential point of 
the Scutum arm. 

te] cite. koyamaKoyama et 
al.(1990) ([cite] cite. koyamaKoyama et al.(1990)) re- 
ported many X-ray transient sources in the Scutum arm 
with a flare of the order of lO'^^ erg s~^, and suggested 
that most are Be/X-ray binaries. Since the hard X-ray 
spectrum, neutral iron line, erratic flux variability and co- 
herent pulsations are characteristic of Be/X-ray binaries, 
AX J1841. 0—0536 would be a member of this sub-class of 
X-ray pulsars. No optical counterpart has been found to 
confirm this classification. The luminosity is consistent 
with that of a "Type I outburst" of Be/X-ray binaries 
(^0 [cite] cite. negueruelaNegueruela(1998)). 

From a semi-empirical relation of the orbital 
period vs. pulse period and the X-ray lumi- 
nosity (§00 [citc]cite.corbetl984Corbet(1984); §00 
[cite]cite.corbetl999Corbet et al.(1999)), we suggest that 
the orbital period of AX J1841. 0-0536 is ~25 days. 
Be/X-ray binaries often exhibit fiares with an in- 
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Table 1. Best-fit parameters for AX J1841. 0—0536 for an absorbed power-law and a Gaussian model.* 

Parameters obs. 1 obs. 2a obs. 2b obs. 2c 

Power-law 

F 1.0 (0.3-1.9) 2.0 (0.9-3.6) 0.7 (0.1-1.3) 1.1 (0.9-1.4) 

Gaussian 

Center (keV) 6.4 (fixed) 6.4 (fixed) 6.4 (fixed) 6.4 (6.3-6.5) 

EWt(xl02eV) (< 4.7) not determined 1.9 (< 4.3) 2.3(1.0-3.8) 

7V| (xlO^^cm-^). 3.2 (1.1-6.4) 3.2 (0.8-8.4) 4.3 (2.0-7.9) 7.2 (6.0-8.6) 

Fluxii (erg cm-^s-^) 2.0 x 10"" 1.9 x lO'^^ 1.3 x 10"" 9.5 x 10"" 

Luminosity§" (erg s^^) .... 2.3 x lO^^ 2.3 x lO^^ 1.5 x lO^^ 1.1 x lO^^ 

xVd.o.f^ 13.2/21 15.1/13 41.9/34 43.2/50 



* Errors and upper limits are at 90% confidence for one relevant parameter, 
t Equivalent width. 

Assuming the solar abundance ratio (^0 [citeJcite.morrisonMorrison, McCammon(1983)). 
§ Estimated in the 2.0-10.0 keV band. 
^ Assumed distance is 10 kpc (see section ^) 

# Degree of freedom. 



terval of the orbital period (c.f. 4U 1223-62, HOO 
[cite] cite. watsonWatson et al.(1982); V 0332-1-53, gOO 
[cite] cite. stellaStella et al.(1985)) due to a gas surge 
during passage near to the peri-astron or across the 
equatorial plane. From AX J1841. 0—0536, we observed 
two subsequent flares separated by ~ 0.6 d, which 
is much shorter than not only the predicted orbital 
period, but also the shortest orbital period known for a 
Be/X-ray binary pulsar, 16.7 days for A0538— 66 (^0 
[cite]cite.corbetl997Corbet et al.(1997)). 

t'eJcite.parmarParmar et 
al.(1989) ([citeJcite.parmarParmar et al.(1989)) re- 
ported that quasi-periodic fiares with a 3.96 hr interval in 
the limited episode of a Be/X-ray binary EXO 2030-1-375, 
which has spin and orbital periods of 42 s and 46.0 day 
respectively (^0 [cite]cite.stollbcrgStollberg(1997)), are 
both different from the flare interval. They inferred that 
the quasi-periodic rapid flares, fast rise and exponential 
decay, are due to instabilities in the either magneto- 
sphere (e.g., ^0 [cite]citc.lambLamb et al.(1977)) or the 
accretion disk (e.g., ^0 [cite]cite.taaml984Taam, 
Lin (1984)). Asymmetric flow of a stellar wind 
(e.g., goo [cite]cite.taaml988Taam, Fryxefl (1988)), 
or non-homogeneity of the OB star wind (e.g., OOO 
[cite]cite.baadeBaade(1984)) may also be a possible sce- 
nario. The same process can be applied to the two flares 
of AX J1841. 0—0536, although the detailed structure, 
the fiare profile, number, size and so on are significantly 
different from those of EXO 2030-^375. 

A remarkable feature of AX J1841. 0—0536 is an erup- 
tive on-set of a fiare with a rise-time of ~1 hour. We 
note that other transient sources, AX J1845. 0—0433 
(^0 [cite]cite.yamauchil995Yamauchi et al.(1995)) and 
XTE J1739-302 (|00 [cite]cite.smithSmith et al.(1998); 
^0 [citeJcite.sakanoSakano et al.(2001)), exhibit multiple 
flares with a fast-rise of a few hours, although no coherent 
pulsation has been found. We suspect that these new tran- 
sients comprise a new sub-class of X-ray binaries, which 
provides a crucial test for the accretion-fiow instability. 



Therefore, further monitor observations on these unusual 
X-ray binaries are encouraged. 
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Fig. 4. Background-subtracted spectra (crosses) and the 
best-fit models (solid histograms). Data residuals are shown 
in the lower panel of each figure. 



